Simulation of noise-assisted transport via optical cavity networks 
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Recently, the presence of noise has been found to play a key role in assisting the transport of 
energy and information in complex quantum networks and even in biomolecular systems. Here 
we propose an experimentally realizable optical network scheme for the demonstration of the ba- 
sic mechanisms underlying noise- assisted transport. The proposed system consists of a network 
of coupled quantum optical cavities, injected with a single photon, whose transmission efficiency 
can be measured. Introducing dephasing in the photon path this system exhibits a characteristic 
enhancement of the transport efficiency that can be observed with presently available technology. 

PACS numbers: 42.50.Ex 03.65.Yz 42.79.Gn 



I. INTRODUCTION 

The presence of noise in quantum transmission net- 
works is generally considered to be deleterious for the 
efficient transfer of energy or classical/quantum infor- 
mation encoded in quantum states. Quantum networks, 
used for the transmission, are unavoidably interacting 
with an external noisy environment and this interaction 
significantly affects the quantum coherence of the sys- 
tem evolution. It is indeed commonly accepted that the 
presence of decoherence [l| is responsible for the unde- 
sired and uncontrolled transfer of information from the 
system to the environment, which in turn reduces the 
coherence in quantum systems. However, recently noise 
has been found to play a positive role in creating quan- 
tum coherence and entanglement @, Motivated by 
fascinating experiments showing the presence of quan- 
tum beating in photosynthetic systems 0-Q , subsequent 
theoretical work pointed to the idea that the remarkable 
efficiency of the excitation energy transfer in light har- 
vesting complexes during photosynthesis benefits from 
the presence of environmental noise [7|, [8[. Indeed, the 
intricate interplay between dephasing and quantum co- 
herence and also the entanglement behaviour during the 
noise-assisted transport dynamics have been elucidated 
in more details in Refs. |9l4l2j. Perhaps even more sur- 
prisingly, the dephasing was recently found to assist the 
transfer of classical and quantum information in commu- 
nication complex quantum networks [13| . 

Recently, quantum optical systems have been exploited 
as a promising platform to simulate quantum processes 
jl34l6|. For example, several implementations of sys- 
tems simulating quantum random walks have been re- 
ported with linear optical resonators 1171 fl8L linear op- 
tical elements [19L [2C|, fiber network [21[ and optical 
waveguides [22|-[25j. Motivated by these results, here we 
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FIG. 1. Analogy between a network of coupled optical cav- 
ities and a physical network where the excitation is carried 
by electrons, such as light-harvesting complexes, (a) The sin- 
gle site of the network is analogous to a single optical cavity, 
(b) The electronic excitation is represented by the presence 
of a single photon in the corresponding optical cavity, (c) 
The transfer of the excitation between two interacting sites 
is analogous to the transfer of the single-photon between two 
adjacent coupled cavities. 



propose a "quantum optical scheme" to investigate the 
noise-assisted excitation transfer process through a set 
of coupled optical cavities. We discuss a four-site optical 
network and derive the set of relevant parameters that 
rule the time evolution of the system. A detailed nu- 
merical simulation of this dynamics, when one cavity is 
injected by a single-photon, is performed employing real- 
istic experimental parameters, showing that the presence 
of a suitable dephasing process in each site of the net- 
work allows for a characteristic increase of the excitation 
transfer efficiency. Furthermore, we consider several as- 
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pects such as phase-stabilization of the cavities and the 
implementation of the dephasing that are necessary to 
observe a clear enhancement of the photon transfer rate 
from one cavity to an external detector mimicking the so- 
called reaction center of the light-harvesting complexes. 
Finally, we investigate how entanglement degrades dur- 
ing the time evolution of the optical network. 

The paper is organized as follows. In Sec. [II] we define 



the model that describes the dynamics of the four-site op- 
tical network analysed in this paper, including the mas- 
ter equation for the two relevant noise processes. Then 
in Sec. IHII we perform a detailed derivation of a realistic 
set of parameters for the system. In Sec. [IVlwe report 
the results of a numerical simulation of the dynamics of 
the network. Finally, the conclusions and final remarks 
are presented in Sec. El 




FIG. 2. Quantum simulation of noiseassisted excitation transfer through a network of coupled optical cavities, (a) Simplified 
scheme of a four-site fully-connected network. The excitation is injected in site 1, and exits from the network by coupling of 
site 2 with the output sink, (b) Equivalent network of four coupled cavities. The excitation is given by a single-photon pulse 
injected in cavity 1. The right output mode of cavity 2 is the output sink channel for the excitation. 



II. MODEL OF THE NETWORK 

In this section we describe in details the model under- 
lying the dynamics of the proposed network of optical 
cavities - a schematic view of this system in relation to 
the light-harvesting complexes is shown in FigJU Start- 
ing from the Hamiltonian describing non-interacting cav- 
ities, one has 

i-Lcav = ^2 fi^&lai , (1) 



with hi and a] being the usual bosonic field operators, 
which annihilate and create a photon in cavity z, and 
uj being the resonance frequency, which we assume for 
simplicity to be equal for all cavities. The transfer of 
photons between the optical cavities is described by the 



following Hamiltonian term: 

where the sum on extends over all the connected 

cavities, and g^ are the coupling constants. Moreover, we 
assume that this dynamics is subject simultaneously to 
two distinct noise processes acting on each optical cavity 
v. 

• a dissipation process that leads to photon loss with 
rate T^, described by the following Lindblad super- 
operator 

£diss(p) = ^2 r i{~ {S-a^pj +2aipaf) , (3) 

i 

• a pure dephasing process that randomizes the pho- 
ton phase with rate 7^ given by a Lindbladian term 
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which has the form 
£de P h(p) = 5Z 7i ( - + 2a\aipa\ai) . (4) 



Besides, the total transfer of the single photon is mea- 
sured in terms of photons detected on the right-hand side 
output of the cavity 2, which represents the so-called sink 
or reaction center of the biological systems (here, a single- 
photon detector), described by the Lindblad operator 

£>Det(p) = T Det [2d ] Det d 2 pd\d Det - {a\a Det a* Det a 2 , p}} (5) 

where a* Det describes the effective photon creation oper- 
ator in the detector and T]j et is the rate at which the 
photon irreversibly gets the detector on the right side 
of the optical cavity 2 - see Figs. [2] and HJ Hence, the 
photon transfer efficiency is measured by the quantity 

Psink 

(t) = 2T Det / Tr\p(t')ala 2 ]dt' . (6) 
Jo 

In the following numerical simulation, we will assume 
that there is a single photon initially in the cavity 1. 
Notice that, since our scheme does not involve any non- 
linear process, a single photon experiment repeated many 
times exhibits the same statistics obtained with an in- 
jected coherent state (26| . 

III. PARAMETERS OF THE NETWORK 

In this section we discuss the experimental details of 
the optical cavity network setup sketched in Fig. [2j in 
order to simulate the mechanisms underlying the noise- 
assisted transport phenomena. The excitation of the net- 
work, i.e. a single photon at wavelength A = 800 nm, is 
generated through a heralded single photon source, based 
on the spontaneous parametric down-conversion process. 
We consider the case of a network of dk = 1 cm long cav- 
ities. The distance between each cavity and the central 
beamsplitter (BS), chosen with transmittivity rj = 0.5, is 
taken to be Ik = 20 cm. More specifically, cavity 1 is cho- 
sen with mirror reflectivities R 1 ™ = 0.9 (for the internal 
mirror pointing towards the BS) and R° ut = 0.99 (for the 
external mirror), cavity 2 with R l 2 n — R™ 1 = 0.9, while 
cavities 3 and 4 with Rf = 0.9 and R° ut = 0.999. The 
loss parameter of each cavity j is given by 

with C6j =0.35 m _1 , while the average number of round 
trips for a photon in the cavity is given by rrij — (1 — 
The parameters adopted for the numerical sim- 
ulation are respectively the set of coupling coefficients 
between the cavities of the network, the dissipation rate 
and the transmission rate from site 2 to the output mode 
k^et (i.e., the detector). 



A. Dissipation rate 

The dissipation rate in each cavity j can be evaluated 
by considering the amount of losses in rrij round trips, i.e. 
the average flight time of the photon in the cavity. Such 
parameter can be evaluated according to the following 
expression: 

where £| represents the fraction of optical power which 
remains confined in the cavity after each round trip, tj 
is the photon flight time in one round trip, rrij is the 
average number of round trips for the photon in the cav- 
ity j, and D is the dissipation in one round trip only 
due to diffraction or coupling with other optical modes, 
i.e. D = 1 — e~ 2oijdj . Eq. (|SJ) evaluates the fraction of 
optical power lost in rrij round trips, divided by the av- 
erage flight time t cav = rrijtj = 2drrij/c with c being the 
speed of light. In our setup, one has Tj ~ bOMHz and 
D ~ 0.007. We consider negligible the losses between 
the cavities and the beamplitter due to the adoption of 
anti-reflection coating optics. Notice that the average 
flight time of the photon in the cavity, i.e. t cav , is of 
the order of ns and defines the time scale of the process. 
Hence, the corresponding linewidth of the cavity alone, 
evaluated from the spectral properties of the intracav- 
ity field, is of ~ 2 GHz. The linewidth of the injected 
photon must be much smaller than this value, hence nar- 
row band parametric down conversion source, such as the 
one obtained through periodically-poled nonlinear crys- 
tals, is necessary for an efficient cavity-photon coupling. 
Let us note that the presence of external mirrors with 
reflectivities i?^* 4 < 1 introduces additive channels T out 
for losing the photon out of the network in spatial modes 
1*1*2*4 ~~ see Fig. El The dissipation rate due to such pro- 
cess is given by the fraction of optical power lost through 
the external mirror in time t cav = rrijtj, and can be eval- 
uated as: 

/i ]Dout\ m i i ( -nout\mn-\-l 

^-^E^'^ — • (9) 

m j z j i=0 m j l j 

For cavities 3 and 4, such additive dissipation channel 
is of the order of T^ t ~ 10 MHz, thus being negligible 
with respect to the dissipation due to intracavity losses. 
For cavity 1, we can introduce a feedback system to dis- 
card those events which correspond to losing the photon 
through this channel. This system is shown in Fig. |4]and 
exploits the polarization state of the photon by inserting 
a Faraday rotator (FR). More specifically, the photon 
with \H) polarization state, after the polarization beam- 
splitter (PBS), is rotated by the A/2 waveplate and by 
the Faraday rotator in the polarization state \V). When 
the photon exits the network from the external mirror of 
site 1, the propagation through the A/2 waveplate and 
the Faraday rotator in the opposite direction maintains 
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FIG. 3. Scheme of all coupling, dissipation and dephasing 
processes of the network. All sites are coupled with the net- 
work with the gkj couplings, and undergo both a dephasing 
process, with rate T, and internal losses, with rate Yj. Site 
1: external losses are reduced with a feedback system. Site 2: 
external coupling with the sink. Sites 3,4: external losses are 
negligible with respect to internal losses Y 3 out <C Yj . 



its polarization state \V) unaltered. Finally, the photon 
is reflected by the polarization beam-splitter and then 
detected. This allows to discard those events when the 
detector clicks, thus allowing to reduce the effective dissi- 
pation term Tl ut . We notice that recent papers reported 
the realization of high detection efficiency (~ 75%) silicon 
avalanche photodiodes and superconducting transition 
edge detectors, with the perspective of reaching a value 
of rjdet ~ 90% [29|, [53. The adoption of these devices 
would reduce the effective dissipation T^ nt due to the ex- 
ternal mirror by a factor 0.2, and hence in our case from 
- 100MHz without post-selection to Yl ut ~ 20MHz. 



coefficients is derived in Ref. [27j and reduces here to the 
expression: 



1 / \ V 2 

1 / V k Vj \ 



(11) 



where vj are the intracavity group velocities, dj are the 
cavity lengths and t^j is the amplitude transmission coef- 
ficient between the fields Aj and A^ in the two sites. The 
amplitude transmission coefficient can be directly evalu- 
ated by analyzing the fields in the classical limit. The 
calculation of this parameter can be divided in three in- 
termediate steps. Particularly, in the following we specify 
these calculations for our setup in Fig. [2j 

1. Output field from cavity k. As a first step, we evalu- 
ate the ratio between the intra-cavity field and the 
field at the output of cavity k and it is given by: 



A 



(k) 

out 



i(*0 



i8cf) k — a k d k 



(12) 



where 5 fa is the phase term due to propagation 
inside the cavity k. 

2. Inter- cavity field at the input face of cavity j. The 
field at the input face of cavity j can be evalu- 
ated as the coherent superposition of all possible 

(k) 

paths of the output field from site fc, i.e. A y Q J t . 
Such quantity strongly depends on interference ef- 
fects between all possible paths that the photons 
can take in the network. We restrict our treatment 
only to the first order path. The ratio between the 
intercavity field and field at the output of cavity /c, 

i.e. h Ai 



int k Y , has the following form: 



B. Coupling rate 

The cross-coupling coefficients have been evaluated fol- 
lowing the theory of Marcuse [27], HH • The time evolution 
of the intra-cavity field amplitude is described by the fol- 
lowing set of first order differential equations: 



.(ny-n)^ + ,5>r^j 



k^j o- 



VG ACT 

Ijk^k > 



V2 



e <4>k+4>i) Kk 



(13) 



dt 

(10) 

where j is the cavity index, (V, a) are the mode indexes, 
s v ° and g%J are respectively the self-coupling and cross- Kkj 
coupling coefficients, is the optical frequency corre- 
sponding to the eigenvalue of the propagation equation 
for the optical mode v in cavity j, Q is a reference fre- 
quency corresponding to the center of the spectrum for 
Qj, and Aj are the field amplitudes. In our case, we con- 
sider only a single-mode of the electromagnetic field and 
the mode indexes (z/, a) can be neglected. An approxi- 
mate expression for the evaluation of the cross-coupling 



where n r is the number of times the photon has 
been reflected by the beam-splitter, fa = 2ii 1j ^i>l is 
the phase-shift due to spatial propagation between 
the cavity k and the BS with vl being the field 
optical frequency and being the distance between 
the cavity k and the central BS, and the form of Kkj 
is different whether the cavities k and j are directly 
linked by the beam-splitter or not, i.e. 

for direct link 
for indirect link 
(14) 

where p and q are cavity indexes satisfying k ^ j ^ 
p ^ q. In order to generalize the expressions of lkj 
and Kkj to the case of a BS with transmittivity 
r] ^ 0.5, each factor has to be replaced with ^/rj 

or y/1 — n whether the photon has been reflected or 
transmitted by the beam-splitter. 
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Intra- cavity field in cavity j. The intra-cavity field 
in site j is related to the intra-network field at its 
input face according to: 



J-'cav 
intra 



(15) 



in the coupling coefficients (induced by the phase fluc- 
tuations) in between experiments, we analyze also the 
case in which there is a static disorder of ~ 20% in the 
coupling rates. 



C. Transmission rate 



hence, the coupling coefficients can be finally writ- 
ten as: 



9k i = 



2i 



dkdj 



1/2 ^(1 -#«)(! 



Rf 



-Oikdk 



(16) 



Following these calculations, the absolute values of the 
coupling rates are found to be g\2 = 4.3 GHz, gis = 5.7 
GHz, #i4 = 7.6 GHz, g 23 = 6.1 GHz, # 24 = 4.5 GHz and 
#34 = 5.9 GHz. To take into account the fluctuations 



Finally, the transmission rate from cavity 2 to the out- 
put mode k on £ is evaluated as above for Tj, by consider- 
ing the amount of field which is transmitted through the 
mirror R™ 1 in the output mode in Vfi2 round trips. In 
other words, the rate at which the photon is transferred 
to the detector can be evaluated with the same expres- 
sion for the external mirror dissipation rate of Eq.(|9]). 
The numerical evaluation of this parameter gives: 



1 



Det 



m^ldjc 



1 GHz 



(17) 




FIG. 4. Experimental setup for the four-cavities optical network. The single-photon in the input mode k; n is injected into 
the network in site 1. The successful transfer of the excitation in the network is given by the detection of the single-photon 
on mode kdet- The coupling coefficients between cavities can be changed by varying the transmittivity and the reflectivity of 
the beamsplitter (BS). The feedback system to reduce the dissipation rate T\ ut from site 1 exploits the polarization degree of 
freedom of the photon as described in the text. [Inset (a)] Sketch of an active phase stabilization apparatus. An auxiliary laser 
is injected and extracted in the network by two dichroic mirrors (DM). The measurement on the auxiliary laser is exploited to 
drive a piezoelectric system which stabilize the cavity lenght. [Inset (b)] Introduction of dephasing rate by modulation of the 
index of refraction through an electro-optical crystal. 



D. Experimental tasks 

The two main challenges for the proposed experimental 
realization, reported in Fig. |4j regard the phase stabil- 
ity of the cavities and the implementation of a suitable 



device to introduce the necessary amount of dephasing 
rate. An accurate control on the cavity length is neces- 
sary in order to maintain the cavities at resonance with 
the photon wavelength and to keep the coupling rates 
constant. This phase stabilization can be achieved by an 
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FIG. 5. Site 1 population behaviour and transfer efficiency 
vs. time (in ns) for the optical cavity network for the noise- 
less (continuous line) and noise (dashed line) case. The pho- 
ton transmission is significantly enhanced by the presence of 
dephasing. The error bars are calculated considering a 20% 
static disorder (10 3 samples). 



active feedback system working on an auxiliary laser su- 
perimposed with the single photon with a dichroic mirror, 
and then measured after its passage through the cavity 
[Fig. HI inset (a)]. 

The parameters of the network here adopted corre- 
spond to a low finesse optical cavity. Hence, a length sta- 
bility of the order of few nanometers is necessary, which 
is a requirement fully achievable with the current tech- 
nology. The dephasing rate inside the cavity can be in- 
troduced by acting on the beam path, through the prop- 
agation factor e lkz . This can be done by modifying the 
propagation length z or by varing the wavevector /c, for 
instance by acting on the index of refraction inside the 
cavity n. The phase modulation can then be achieved 
by different methods, depending on the desired modula- 
tion rate. For a dephasing rate of the order of 1 GHz, as 
the one exploited in the numerical simulation below, an 
electro-optic or acusto-optic modulator can be inserted 
in the cavity and exploited to shift the frequency of light 
[Fig. IU inset (b)]. 



IV. RESULTS OF THE NUMERICAL 
SIMULATION 



In this section we investigate the dynamics of the opti- 
cal cavity network, by using our model with the coupling 
rates gij as estimated above, the dissipation rate equal 
to Tj = 70 MHz for all the sites by taking into account 
the different dissipative processes, and the transfer rate 
to the detector being Tjj et = 1 GHz. In particular, in 
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FIG. 6. Dependence of p S ink at a fixed time t = 20 ns as a 
function of the dephasing rate 7. These results show the re- 
markable robustness of this process, supporting the possibility 
of being experimentally observed. 



Fig. [5]we show the cavity photonic population as a func- 
tion of time in both cases of local dephasing (with rate 
1 GHz) and no dephasing. The introduction of dephas- 
ing reduces the destructive interference between all the 
possible photon pathways and increases remarkably the 
overall transfer efficiency from about 40% to more than 
70%. In absence of dephasing, the photon is trapped in 
some superposition (dark) states which are not coupled 
to the site 2 and this explains the lower transfer efficiency 
- see Ref. [lOj for more details. Moreover, we consider 
also the case in which the coupling rates suffer a static 
disorder of 20% and also the dependence of the transfer 
efficiency as a function of the dephasing rate - see Fig. 
[6l These results further supports the fact that the noise- 
assisted transport could be experimentally observed by 
the present optical setup. 

Finally, to quantify the entanglement dynamics, we 
study logarithmic negativity [31(, i.e. E(A\B) = 
log 2 ||p rA ||i, measuring the entanglement across a bipar- 
tition A I B of a composite system, where Ya is the par- 
tial transpose operation of the density operator p with 
respect to the subsystem A and || • ||i denotes the trace 
norm. In Refs. [To|, Hi) , in the context of the light- 
harvesting systems, it was found that the increase in the 
transfer efficiency is not strictly related to the presence 
of entanglement between the sites of the network and 
a similar behaviour has been numerically found in this 
cavity dynamics. However, to explore some more the 
capabilities of this optical cavity network as a conduit 
for not just energy (or classical information) but quan- 
tum information, we show in FigJ3 how another form 
of entanglement (in this context more relevant, as mea- 
sureable more directly) degrades through it. To that 
end, we introduce an ancillary photon, which initially 
shares a maximally entangled state with the single pho- 
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FIG. 7. Degradation of entanglement between two photons 
(initially in a maximally entangled state), in which only one 
photon is sent to the network through the cavity 1. It is mea- 
sured in terms of log- negativity between the ancilla photon 
and the photon leaving from the cavity 2 (no detector). 



V. CONCLUSIONS AND OUTLOOK 

We proposed a quantum optical network, based on a 
set of coupled cavities, in order to investigate the effects 
of noise in the excitation transfer. A detailed numeri- 
cal simulation for experimentally realistic values shows 
the presence of a substantial enhancement in the photon 
transport efficiency when dephasing noise is introduced 
in the cavity. As a final remark, we note that a similar 
scheme can also be implemented by exploiting a network 
of atoms in suitable ion traps [34]. The results reported 
here may open interesting perspectives for a deeper inves- 
tigation of the fundamental mechanisms that underly the 
very high efficient excitation transfer in light harvesting 
complexes and for possible applications in solar cell tech- 
nology. Finally, this type of experiments can also trig- 
ger significant activity on two different areas, namely the 
modeling of complex environments via controlled inter- 
actions and the development of noise-assisted protocols 
for quantum communication [l3j |. 



ton injected in the optical cavity 1, i.e. in the EPR state 
1/V2(|0)anc|0)i + |l) onc |l)i), with |0) and |1) represent- 
ing, respectively, the absence and presence of a photon 
[32l . [33j . As the system evolves, the entanglement be- 
tween the ancillary photon and the photon leaving out 
from the cavity 2 oscillates in time and is almost vanish- 
ing in presence of dephasing. 
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